Stationary Motion of the Adiabatic Piston 
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We consider a one-dimensional system consisting of two infinite ideal fluids, with equal pressures 
^"*^ ■ but different temperatures T\ and Ti, separated by an adiabatic movable piston whose mass M is 

much larger than the mass m of the fluid particules. This is the infinite version of the controversial 
adiabatic piston problem. The stationary non-equilibrium solution of the Boltzmann equat ion for 
the velocity distribution of the piston is expressed in powers of the small parameter e = yJm/M, 
and explicitly given up to order e 2 . In particular it implies that although the pressures are equal 
• on both sides of the piston, the temperature difference induces a non-zero average velocity of the 

piston in the direction of the higher temperature region. It thus shows that the asymmetry of the 
fluctuations induces a macroscopic motion despite the absence of any macroscopic force. This same 
conclusion was previously obtained for the non-physical situation where M = m. 



I. INTRODUCTION 
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The adiabatic piston problem has been the subject of continuous controversy for many years (see [1] and references 
given therein). The problem is the following. The system is a finite cylinder containing two gases separated by an 
adiabatic movable piston. Initially the piston is fixed by a brake and the two gases are in equilibrium. At a certain 
time the brake is released and the question is to find the final equilibrium state. Formulated differently, it is one of 
the "Problems in statistical mechanics that I would like to see solved", mentionned by E. Licb in his talk at Statphys 
20 [2]. The controversy is related to the fact that it is a well known example, discussed in many text books (e.g. [3]), 
for which the principle of entropy maximum, i.c the 2 nd law of thermostatics, is not sufficient to find the equilibrium 
iy-j ' position. Indeed the equilibrium position depends in an essential manner on the viscous force and one has to solve 
explicitly the time evolution to find the final position (just like in mechanics when solid friction is present). In other 
words it is a problem which can be solved within the framework of thermodynamics (of irreversible processes) , but 
not within the framework of thermostatics [1] . 

In this context the question has often been raised whether the adiabatic condition is an ideal constraint, which is 
essential to discuss the foundation of thermodynamics [4] , but which can never be realized in a real physical situation 
(just as the condition of a frictionless piston). The question has even been extended to whether the 2 nd law does not 
apply in certain situation. 

As we know, heat conduction through a wall is associated with the internal motion of the molecules of the wall. If 
the wall is a rigid ideal solid which is fixed, the atoms of the gas will make purely elastic collisions on the wall and 
it is possible to maintain at different temperatures two gases separated by such a wall. The wall is thus adiabatic. 
The problem is then to analyse what happens when this same ideal wall is no longer fixed, but is allowed to move. 
The question is whether the stochastic motion of the wall, induced by the fluctuations in the force exerted by the 
particules of the gas, will lead to an equilibrium state in which the temperatures of both gases are equals. In other 
words whether an adiabatic fixed wall becomes diathermic when it is allowed to move. 

Recently this problem was investigated using the following model. The adiabatic piston was considered as a heavy, 
structureless, particule of mass M submitted on both sides to purely elastic collisons with the atoms of two gases 
in equilibrium with temperatures T\ and T2 respectively [5]. It was assumed that both gases are ideal and made of 
identical point particules of mass to. Considering first that M 3> to, a linear Boltzmann equation was derived for the 
velocity distribution of the piston. The rigorous stationary solution was then found for the special case where one 
takes M = to in the Boltzmann equation and the two gases extend to infinity, respectively to the left and to the right 
of the piston. The main conclusion was that the fluctuations conspire with the spatial anisotropy (Ti ^ T2) to induce 
a macroscopic motion of the piston, with finite average velocity in the direction of the higher temperature region, 
although no macroscopic force is present (the pressures on both sides of the piston are equal). The same model with 
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M = m was then considered in [6] for a system with finite length; it was shown that the piston evolves toward an 
equilibrium position corresponding to a uniform density of particules throughout the system. 

In the following we study the more physical situation where M 3> m. In this case the Boltzmann equation describes 
asymptotically exactly the dynamic of the piston. This is because the massive piston has a vanishing probability to 
interact back with the perturbation it causes in the states of the surrounding fluids (no recollisions in the sense of 
the kinetic theory) . These perturbations go away to infinity, and the piston at consecutive collisions always " sees" on 
both sides the unperturbed initial equilibrium states. This same model was alredy considered by J.L. Lebowitz [7]. 

In the special case of equal temperatures on both sides of the piston, i.e T± = T2, the problem reduces to the 
classical Brownian motion of the "Rayleigh piston", which has been extensively studied in the literature [8-11], 
using for example the Kramers-Moyal expansion in powers of 1/vM [12]. The general problem of a heavy particule 
undergoing elastic collisons with an infinite ideal gas of light atoms in equilibrium has been rigorously treated by 
Holley in one dimension [13], and by Diirr, Goldstein, Lebowitz [13], followed by Goldstein, Guetti [14] in arbitrary 
dimension. In the limit where m — > and the density of the gas in proportional to mT 1 / 2 , these authors established 
the convergence of the mechanical process describing the velocity of the heavy molecule to the Ornstein-Uehlenbeck 
process. Therefore when the temperature of the bath is uniform the heavy particule will finally thermalized, acquiring 
the temperature of the surrounding thermostat. The characteristic velocities of the piston and of the gas particles 
would be yJksT/M and yJk^Tjm, respectively. So, as a natural small parameter for our investigation with T\ 7^ T 2 
we choose 




The mechanical model is presented in section II. In section III we develop the perturbative analysis of the stationary 
solution of the Boltzmann equation for e — > 0. Section IV contains the evaluation of the average velocity of the piston 
and concluding comments. 



II. MECHANICAL MODEL 



We consider a gas composed of identical point particles of mass m moving in R 1 . The motion between collisions is 
free. When two particles collide they instantaneously exchange their velocities, the collision being perfectly elastic. 

Up to the initial moment t = t a particle of mass M, fixed at some point, isolates the volume of the gas to its 
left (subsystem 1) from that to its right (subsystem 2). Collisions with the fixed mass M are assumed to be elastic 
(specular reflections). The two infinite volumes of the gas on both sides of the mass M are at thermal equilibrium 
characterized by the number densities n±, rii and temperatures Ti, T%, respectively. The equilibrium pressures are 
given by the perfect gas equation of state 

Pi = nik B T u i = 1,2, (2) 

where ks is the Boltzmann constant. 

At the moment t = to the constraint fixing the mass M is removed, and the released mass starts to move freely 
suffering elastic collisions with the gas particles. If before the collision the masses M and m have velocities V and v, 
respectively, their velocities after the collision are given by 

m + M 

v >= v +——(V-v) 
m + M 

Our object here is to determine the asymptotic stationary velocity distribution of the mass M in the case where 
M » m, assuming the initial pressures (||) to be equal 

Pi = P2 = P (4) 

From the physical point of view we thus consider an infinite volume version of the adiabatic piston problem discussed 
in the introduction and which is essentially a one-dimensional problem. The massive structureless particle plays here 
the role of the movable adiabatic wall (or piston) , and the fundamental question concerns the nature of the asymptotic 
stationary state. The assumed equality of the initial pressures (^|) represents the compensation of the momentum 
fluxes arriving at the fixed piston 
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pi = I dv v 2rnv4>™ 1 (v) = / dv v (-2mv)cb™ 2 (v) = p 2 (5) 

JO J -oo 

where <$p denotes the Maxwell distribution 

The average motion of the released piston will thus reflect the effects of fluctuations in the momentum transfers at 
binary encounters. 



III. PERTURB ATIVE ANALYSIS OF THE STATIONARY STATE 



As discussed in [5], the velocity distribution function <$> t (V) of the piston in the stationary state with equal pressures 
on both sides is given by the solution of the linear Boltzmann equation. With the use of the Dirac ^-distribution we 
can write it in the form 



m + M 



(v 1 -V')-V 



dv' J dv'(v' - v')d(v' - V)y4>tiW)MV')s (v 

- J dv(v - V)6(v - V)±-<t% (v)* c (V) 
dV J dv'(V - v')0(V - J)jrW a (v')MV')S (V + -^jjW V) - v\ 



(7) 



dv(V-v)6(V-v)—^ 2 (v)^ e (V) = 

J-2 

Equation (]?]) describes the situation where the temperatures T\ and T 2 characterize the initial states to the left and 
to the right of the piston, respecti vely. 

With the small parameter e = \J m /M introduced in section I we have from equations (|^) 
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So, the integration over V in (0) yields the equation 
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dv (v-V)-4%M^ V ) 



dv(V-v)-<f>%(v)$ e 

±2 



dv(V - v)-^ 2 (v)MV) 

J-2 
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(v-V) 
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V+- s{V-v) 



1-e 2 
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(9) 



which can be conveniently rewritten in the form 
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1-e 2 



l + e' 
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dv(v - V)-mv) + / dv(V - v)-<f>%(v) MV) 



To 



1-e 2 



l+_£ 

1-e 2 



2\ 2 r V 



dv(v - V)—<t%(v) {$ e [V- ^— 5(1; - V)] - MV) 



( V -V)]-* e (V) 



dv(V-v) ¥ ^ 2 (v)\^[V 
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2e 2 
1-e- 



(10) 



Using the expansion of $ e [V — 2e 2 (v — V)/(l — e 2 )] in powers 
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we obtain 



2e 2 (v-V) 



= 2 



(l + e 2 ) 2 E 



(l + e 2 ) 2 E 



dv(v-V)±-<f>% 1 (v) + J dv(V-v)^ 2 (v)\^ e (V) 

^ (-1)" (2 £ 2 ) 
^ n! (1 - e 2 )" 

n=l v ' 

2\n-l i-V 



J — oo 
n (r, 2\n— 1 /-oo 



n+l 



1 / A \ n 



(V) 



1 (2e 2 ) 



^ n! (1 - e 

n— 1 x 
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The series in (|12|) corresponds to the e-expansion of the collision law (|3|). 
Introducing the functions I n (V) defined by 



dv(V-v) n —^(v)-(-l) n / dv(v-VT-cjf T \{v) 



T 2 



we get from (h_2|) the following equation for $ e (V) 



= 2I 1 (V)MV) + (l + e 2 ) 2 ^^^(n + l)I n+1 (V)^J Z e (V) 
Let us first establish some useful properties of the functions I n (V). From the definition we 
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We define Z_i by 



- /n-xCV), n = 0,l,2,. 
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mV 1 



exp 



mV 2 
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It will be convenient to introduce the expansion of I n (V) in powers of V 
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Using equation (15) yields 



In .q — In—l,q—l ^ 
I n ,q = In-q,Q, Vq < n 



I, 



I-l.q-n-1, Vg > n + 1 



From (|16|) and (19) we obtain for q > n + 1 



ln,n+2l 



and 



(-1)' / m ^ i+1/2 



^n,n+2;+l — 7=- 777 
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On the other hand, for q < n, the relation (19) yields I nq = I n - q (V = 0) and 



(1/ \ (™-<?)/ 2 



2""9+ l r(l + (n-q)/2) 



with 
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T(l + l) = l\ 
2k + l _ _^/tF 
2 ' ~~ 2 k ~ 



(23) 



(24) 
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Using the relation (iq), one can rewrite equation (14) in the form 



d j / 2e 2 



(2 + n + ne^ 2+ „(U) ( — ) $ e (^) > = 



(25) 



Integrating (|25|) over the interval [V, oo), and assuming that the integrand vanishes in the limit V — > oo, we obtain 
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(2 + n + ne^)7 2+ „(F) ( — ) $ e (F) = 
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For Ti = T 2 = T, the solution of (0), or @, is the Maxwell distribution 

^e(V)\ Tl=T2 =T = <${V) = 

This suggests to introduce the dimensionless variable 
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and the normalized distribution F t (x) defined by 

F e (x) -- 



AI 



AI 



(26) 



(27) 
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We thus have 



(30) 



Using the expansion ([17]) for l2+ n (V) and the change of variables ( |28| ) we find that the distribution F e (x) is the 
solution of the equation 
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If (n — g) > —2, equation (|23| ) implies 
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-n— 3 



M r(i + 2±jp) 



while for g = n + 3 + 2/, Z = 0, 1, 2, we have from (E2n 



(31) 



(32) 
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We thus obtain the following equation for F e (x) 

oo oo n-\-q q / j \ n 

EE7T^yr( 2 + - + - 2 ) 2 " J «-^ U Fe(a;) = (34) 

where in terms of the parameter 

/ rp x i/4 

|1 (35) 



we have 



\r r - (-t)- t '] 

J r = j± -. — \ ' { , r , for r > -2 (36) 
2 r+3 r(l + (Y + 2)/2) ~ 

J r = [r r + r- r ]^-^-, for r = -3 - 21, 1 = 0,1,2,... 
v 7r 

J r = 0, for r = -4 -22, / = 0,1,2,... 

Let us remark that since Jo = 0, the left hand side of ([mJ) vanishes for e = 0, which is consistent with the assumption 
made when deriving equation (^6|). Therefore, we can divide ( fji] ) by e. Finally, assuming that F e (x) has an asymptotic 
expansion 

oo 

F e (x) = J2z l Fi(x) (37) 
;=o 

in a series of integrable functions Fi(x), we are led to solve at successive orders in e the equation 

°° °° °° °° (<n - 1 +tnll t-9 f d\ n 

EEEE^ M -'^(H«»')n,-i i HW _ (38) 

n=0 g=0 m=0 2=0 v ' * \ / 

The normalization in (|3^) is supposed to be entirely contained in the zero order term 



dxF e (x) = / dxF (x) = 1 (39) 
We will thus require that 

for r = l,2,... (40, 

To obtain the solution for ([38]) at successive orders r = 0, 1, 2, we give below the explicit form of J r 

Ji = {t + t-^/U^, J- 1 = (t + t- 1 )/2V^ (41) 
J 2 = (r 2 + r- 2 )/26, J_ 2 = (r- 2 - r 2 )/2 
•h = (t 3 + r- 3 )/12(V¥, J_ 3 = (r 3 + r- 3 )/V^ 



The terms in (38) which contribute to the order e r correspond to the values of the summation variables n,q,m,l 
satisfying the relation 

n + q + 2m + l-l = r (42) 
The last of equations (E^) shows that some sets satisfying yield a vanishing contribution. 
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Using equation (36) we can now construct the expansion (37) in a systematic way. 
Order e° 
When r = 0, we find 



2 J^x + 6 Ji 



d.T 



F (z) = 2J_! 



l__d 

2 dx 



F (x) = 



Hence the normalized solution reads 

F (x) = exp(-x 2 )/y/ir 
We thus have at the order zero (see equation (gjj)) a gaussian probability density 

MV 2 



M 



2irkB\ / TiT 2 



exp 
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Order e 



When r = 1, equation (Bq) yields a relation between -Fi(a;) and -Fo(^) of the form 



2 J_ix + 6 Ji 



dr 



Fi{x) + 



J- 2 x 2 + 16 J 2 



Using (Pq) we can rewrite (E6J) as 



1 d 



1 / d 



2 4 V ds 



F (x) = 



F (x) 



The condition 



dxFi{x) = 



determines the solution in a unique way. We find Fi(x) = Ai(x)F (x), where 

A 1 (x) = v ^(r-r- 1 )(^-^ 3 ) 
So, up to terms of order e, the distribution F e reads 

F e (a:) 



1 + c\/¥(t - r _1 )(a; - -a; 3 ) 
o 



F (x) 



(43) 
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(48) 



(49) 



(50) 



The general frame developed here permits to continue the determination of higher order terms in a systematic way. 
The calculation of the second order contribution F2(x) is presented in the appendix. 



IV. AVERAGE VELOCITY OF THE PISTON: CONCLUDING COMMENTS 



Using the definition (|29j) we can evaluate the average velocity of the piston in the stationary state associated with 
different temperatures Ti, T 2 , but with the same pressures on both sides. We find 



< V M I dxxF *( x ) 



(51) 
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— eV^(r-r )[< x > -- < a; 4 >] 
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Equation (^l|) contains the fundamental prediction of our analysis. It shows that the average velocity of the piston is 
different from zero although the pressures are the same. Moreover, if T2 > Xi, the velocity is positive showing that 
the piston evolves in the direction of the high temperature region. In conclusion we thus recover in the more physical 
situation of M ^S> m, qualitatively the same result as in the equal mass case M = m studied in [5]. The advantage 
of the present analysis is that the Boltzmann description is asymptotically exact when e = yro/M — * 0, because the 
effect of recollisions vanishes in this limit. We thus conclude that for a massive piston a macroscopic motion is induced 
despite the equality of the pressures, reflecting the asymmetry of fluctuations. This result, obtained for an infinite 
volume, supports the point of view that the adiabatic piston separating finite volumes of gases at equal pressures but 
different temperatures will be also systematically displaced until the uniform distribution is attained (for a rigorous 
analysis of the special case M = m see [6]) The average kinetic energy of the piston is determined by the zero order 
Maxwell distribution (|3),(B|), with temperature \/T{F 2 . We thus have 



M y2 k B ^T{r 2 
2 2 



(52) 



It follows from (|l|) and (47), that 

< V >-- 
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(53) 



Let us end this section by writing down the asymptotic formula for the stationary distribution $ e (V). When e 

mV 3 



s/mfM < 1, 
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6fcsv / 7i 7 Z2 



(54) 
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APPENDIX 

Order e 2 

From the expansion (38) there follows an equation relating F^ix) to Fq{x) and F\{x) 

[2J_ lX + 6J 1 —]F 2 (x) + [J_ 2 x 2 + 16J 2 f — J }F 1 {x) + (55) 

3 \dx J dx \dx J dx 

Putting F 2 (x) = A 2 {x)Fq(x) we get with the help of ( |36| ) a first order equation 

j-A 2 {x) = (r - r- 1 ) 2 [7r(i a; 5 - yx 3 + 4a:) + (yx 3 - te)] (56) 



The solution 



A 2 {x) = (r - r- 1 ) 2 ^ + (2tt - 4)a; 2 - 1(5tt - 7)x 4 + ^-irx 6 } (57) 

6 18 



contains a constant A which can be determined from the condition (ECM). One finds 



93 Q 

< 58) 

Hence, the second order contribution to F e (x) is given by 

93 Q 1 1 

F 2 (x) = (r - ^H-tt - - + (2tt - 4).x 2 - -(5tt - 7)x 4 + -^ 6 ]F (o;) (59) 

Let us remark that the general structure of the series expansion ( p7[ ) has the product form of a polynomial in variable 
x multiplying the gaussian distribution Fo(x) (]44|). At the order ^ one gets terms in x 5 ,x 7 ,x 9 , at the order e 4 terms 
in a; 8 , a; 10 , a; 12 appear. In general, the contribution of the order e r yields a polynomial factor in a; 3r " 4 , a; 3r ~ 2 , x 3r . 
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